Isotopic shift in angle-resolved photoemission spectra of Bi2Sr2CaCu208 due to 

quadratic electron-phonon coupling 
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In connection with the experiment on oxygen isotope effect of Bi2Sr2CaCu2 08 with the angle- 
resolved photoemission spectroscopy (ARPES), we theoretically study the isotope-induced band 
shift in ARPES by the Hartree-Fork and quantum Monte Carlo methods. We find that this band 
shift can be clarified based on a quadratically coupled electron-phonon (e-ph) model. The large 
ratio of band shift versus phonon energy change is connected with the softening effect of phonon, 
and the positive-negative sign change is due to the momentum dependence of the e-ph coupling. 
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I. INTRODUCTION 

The study of superconductivity has a long history. For 
the conventional superconductors, the mechanism has 
been successfully explained by the BCS theory^, that the 
electron-phonon (e-ph) interaction is at the central stage. 
For the high-Tc superconductors, however, the mecha- 
nisms is still not clear with a number of basic questions 
need further investigation. In order to clarify these prob- 
lems, enduring efforts have been made for two decades 
by the experimentalists and theorists. 

In the experimental aspect, since the angle-resolved 
photoemission spectroscopy (ARPES) directly measures 
the electronic occupied states in a crystal, it has be- 
come one of the most important experimental techniques 
to study the electronic properties of superconductors. 
Nowadays, with the great improvement of energy reso- 
lution, ARPES is able to observe the fine structures due 
to e-ph interaction^ and the tiny superconducting gap^ 
with an energy resolution less than 1 meV. Recently, the 
oxygen isotope effect has been studied with ARPES on 
the high-Tc superconductor Bi2Sr2CaCu208 (Bi2212) by 
two groups^. Comparing their data, at least one com- 
mon feature has become clear, the spectra are shifted 
when ^^O is substituted by its isotope ^^O, whose move- 
ment in this crystal is regarded as a lattice vibration, 
i.e., phonon. This isotope induced band shift indicates 
that the e-ph interaction has an effect on the electrons, 
hence providing direct evidence for the interplay between 
electrons and phonons in this material. 

Since the first report by Gweon et alA, the isotopic 
band shift has become an controversial issue and gained 
considerable interest^i^, as the observed band shift is up 
to 40 meV, much larger than the maximum isotopic en- 
ergy change of phonon (^ 5 meV) according to the mea- 
sured vibration energies of oxygen^. In addition to this 
large ration of band shift versus phonon energy change, 
the presence of both positive and negative band shifts in 
the ARPES also implies that the electrons and phonons 
are coupled in a complicated manner. Very recently, Dou- 



glas et al^ repeat the experiment and find the shift is 
only 2±3 meV, which is inconsistent with the large shift 
found by Gweon et alA. 

In this paper, we look into the e-ph interaction in 
Bi2212 and its relation with ARPES by the Hartree-Fork 
and quantum Monte Carlo (QMC)^*' methods. Our pur- 
pose here is to figure out whether the isotope-associated 
anomalies, especially the large ratio of band shift versus 
phonon energy change, are possible at all from a theoreti- 
cal point of view. We will show that the abovcmcntioncd 
band shift can be explained within a scenario of phonon 
softening effect driven by the e-ph coupling, while the 
positive-negative sign change is due to a momentum de- 
pendence of the coupling. The remainder of this paper 
is organized as follows. In Sec. 2, we put forward the 
model Hamiltonian. In Sec. 3, an investigation on this 
model with Hartree-Fork theory is conducted. In Sec. 4, 
we show the numerical results of Hartree-Fork and QMC 
methods, and discuss the origin of band shift. Our con- 
clusion is presented in Sec. 5. 



II. MODEL 

In the theoretical aspect, to investigate the effect of lat- 
tice vibration on the electrons, two kinds of e-ph interac- 
tion are usually modeled for cuprates, i.e., diagonal and 
non-diagonal couplings^ The former is a correlation be- 
tween the phonons and intra-site electron densityi^, and 
the latter is between the phonons and inter-site chem- 
ical bondi^. In the Cu02 plane of cuprates, as shown 
in Fig. 1, when an electron is hopping from one copper 
atom to another, its movement is strongly affected by 
the vibration of oxygen atoms between the initial and fi- 
nal sites^'*. This kind of interaction thus turns out to be 
an off-diagonal type coupling, which modulates the elec- 
tronic band structure. In this sense, the off-diagonal e-ph 
coupling is crucial to the aforementioned isotope-induced 
band shift. 

In this work, since we are primarily concerned with the 
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FIG. 1: A Schematic plot of Cu02 conduction plane in 
cuprates. The copper atoms (black circles) form a simple 
square lattice, and the oxygen atoms (white circles) are lo- 
cated between the nearest-neighboring Cu sites. 



effect of band shift, we start from the following model 
Hamiltonian {h = 1 and A;b = 1 throughout this paper): 

H = -^^^t{l,l + Si){al^ai+s^,a + al_^_s.^^ai,a) 

1,(7 i,Si 

l,<y (1,1') ^ " ^ 

where (a;<j) is the creation (annihilation) operator of 
an electron with spin a at the Cu site I on a square lattice. 
t{l,l + 5i) is the electronic transfer energy between two 
Cu sites I and I + Si, and Si enumerates the ith {i=l or 2 
in this work) nearest neighboring Cu sites, which is acces- 
sible from site L /i is the chemical potential, defining the 
total electron number in this system. The oxygen phonon 
is assumed to be of the Einstein type with a frequency wq 
and a mass m. A (= 1 + Am/m) is the mass change factor 
of phonon due to the isotope substitution. In the third 
term, qw is the dimensionless coordinate operator of the 
oxygen phonon located between the nearest-neighboring 
Cu sites I and I' , and the sum denoted by {I, I') just means 
a summation over all the phonon sites in the lattice. 

In the conduction plane of Cu02, the electronic hop- 
ping integral t{l, I') can be expanded to the second order 
terms with respect to the phonon displacements. For ex- 
ample, in Fig. 1, the nearest-neighbor hopping t{l,2) 
and the second-nearest-neighbor hopping f(l,4) can be 
expanded as, 



t{l,2) = h + siql^, (2) 
t(l,4) = t2 + S2{ql2 + qL + <li3 + 44)- (3) 

Here we note that the linear e-ph coupling does not ap- 
pear in the expansion owing to the lattice symmetry of 
present model. Therefore, if we consider the hopping up 



to the second nearest neighbors, Hamiltonian (1) can be 
rewritten as, 

l,a Si 



I, a- $2 



t2 + S2 ^{qll+S, + <lf+Sul+S2) 

Si 

where ti {t2) is the bare nearest- (second-nearest-) neigh- 
bor hopping energy, and si (52) is the off-diagonal e-ph 
coupling strength due to the nearest- (second-nearest-) 
neighbor hopping effect. The sum denoted by Yl'si 
the second term of Eq. (4) is over the common near- 
est neighbors shared by the Cu sites I and I + S2- Here 
we want to stress that by rising this model, we have as- 
sumed the anti-ferromagnetic ordering is established as 
a background. Accordingly, the double occupancy does 
not occur and the on-site Coulomb repulsion between the 
electrons is omitted. Mean while, the inter-site electron- 
electron interaction is partially included in the screened 
values of t's and s's. 



III. HARTREE-FORK TREATMENT 

In order to have an insight in the effect of quadratic e- 
ph coupling, we first study the model Hamiltonian (4) via 
a Hartree-Fork treatment. By applying the Hartree-Fork 
approximation (HFA), 



AB^A{B) + B{A)-{A){B), 



(5) 



we can separate the e-ph coupled term in Eq. (4) and 
get an effective Hamiltonian as. 



HuFA 



He + -ffph, 

X (4<T«;+'5i,<^ + «i+5i,a«i,<^) 
-^EE[^2+4S2(«2(A))] 



(6) 



Wo \ ^ 

2 ^ 

{1,1') 



(7) 
(8) 
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where (• • •) means the expectation value of an operator, 
and 



27Vu;o I V 

L.<T \ Oi 



(9) 



As win be shown later, 7 acts as a renormalization factor 
for the phonon energy. 

According to the many-body theory, the spectral func- 
tion can be obtained from the single-electron Green's 
functioni^. After some algebra, we get the electronic 
spectral function from Eq. (7) as. 



(10) 



where E\^{X) is the phonon- mediated tight binding energy 
of an electron with momentum k, 

-Bk(A) = -2[ti+ si(q'^ {X))]{cos kx + cos ky) 

-4[t2 + 4s2((7^(A))] cos cosky - fi. (11) 

Hence, the isotope induced band shift Ai?k [= Ei^{Xo) — 
E-k(X)] has the following form, 

A^k = -SkA{q^), (12) 

where, 

Sk — 2si(cos kx + cos ky) + 16s2 cos k^ cos ky, (13) 

is the momentum dependent coupling constant, and 

A(q') ^ (q'iXo)) - (q'W), (14) 

is the phonon spatial variation due to the isotope sub- 
stitution. Here, one can clearly see that the band shift 
Ai?k is highly anisotropic due to the momentum depen- 
dence of coupling Sk. This anisotropy can give rise to 
the above-mentioned sign reverse of the isotope effect. 
We shall return to this point in the next section. 

Since the e-ph coupling is in a quadratic manner, the 
phonon wave function retains a Gaussian form, and we 
assume it to be 



D ( 
-^exp' 



(15) 



where Z? is a coefficient to be determined by the variation 
method. By using the Gaussian integral formula. 



dqeyiY>{—aq 



(16) 



/ 1 / TT 
dqq^ exp(-ag^) = - W — , (17) 
-00 z V Ck 



we have 



9' , . , D'' 



2 ' 



and the phonon energy i?ph (= {(t>w\Hpi,\(f)u>)) is 
Nujo f D"^ 1 - 7' 



E, 



ph 



A L»2 



(18) 
(19) 

(20) 



Minimizing the phonon zero point energy of Eq. (20) 
with respect to Z?, the coefficient D is finally determined 



D=t/{l-l)X. 



(21) 



Substituting Eq. (21) into Eq. (18), the phonon spatial 
variation is found to be 



A(g2) = 



1 



2yr^ 



1 



/Af^ 



1 

71 



(22) 



Thus, within the HFA, we can self-consistently determine 
the isotopic band shift in terms of Eqs. (9), (12), (22). 

In terms of Eqs. (20) and (21), the phonon energy 
becomes 



E, 



Nlo, 



ph 



VA' 



where, 



7, 



(23) 



(24) 



is the effective phonon energy renormalized by 7, which 
means the phonon is softened when the e-ph coupling is 
switched on. This phonon softening effect is not due to 
the photoemission or isotope replacement, but the intrin- 
sic e-ph coupling which is in the quadratic manner. Here 
we want to stress that the bare phonon energy luq is not 
an experimental observable, meanwhile the effective one 
LUs can be detected by the Raman spectroscopy. We now 
introduce the isotopic phonon energy change as 



AlO = LUs{Xo) - UJs{X) 



(25) 
(26) 



From Eqs. (12), (22) and (26), we find the ratio AEi,/Au; 
can be represented as 



AEu 



Sk 



Auj 2a;o(l-7)' 



(27) 



which is only subject to the e-ph coupling. In this sense, 
A£^k/Aa; can be regarded as a quantity indicative of the 
strength of e-ph coupling and the degree of phonon soft- 
ening in the e-ph system. 
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FIG. 2: Energy dispersion of Bi2212 from the tight-binding 
fit of this work (solid curve) and Ref. 16 (dashed curve). 
(Parameter selection is described in the text.) The inset shows 
a quadrant of the Brillouin zone with Fermi surface. 
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IV. RESULTS AND DISCUSSION 

As mentioned above, the most intriguing behavior of 
the isotope effect is associated with the ratio AE'k/Ao;. 
We shall focus on this ratio in the present section. 

In order to compare the HFA results with the ex- 
act calculation, we employ the QMC method^ in this 
work. Starting from Hamiltonian (4), we first calculate 
the spectral function A\^{\ cj), and then determine E-^(\) 
by the moment analysis of the spectral function as. 



£^k(A) = / LjAk(A,cj)dw. 



(28) 



In QMC, the temperature is set at T=0.05. Since this 
simulation suffers from the random error, in the numeri- 
cal calculation we impose a large isotopic mass enhance- 
ment to suppress the fluctuation, i.e., from Ao = l to A=2. 
It should be noted here, although this mass variation is 
larger than the ^^0/^*0 substitution, it actually has no 
effect on the ratio AiJk/Aw, which is mass- independent, 
as shown in Eq. (27). 

In Ref. 16, five tight binding parameters are suggested 
for a best agreement with the experimental band struc- 
ture and Fermi surface of Bi2212. While, in this work, 
since we are interested in the mechanism of band shift, we 
restrict our discussion in a case with only nearest- and 
second-nearest-hopping effects, and thus introduce two 
tight binding parameters, i.e., the bare electronic trans- 
fer energies t\ and ti. To reproduce the band structure 
and Fermi surface, we set t\—Q.Tl eV and make it the 
unit of energy for the numerical calculation. In the unit 
of ^1, we assume i2=-0.2, and the chemical potential [i—- 
1 which is invariant with the isotope substitution. The 
bare phonon energy is set as wo=l-Oi and the two e-ph 
coupling constants are fixed at a ratio of si:s2=l:-l, for 
simplicity. In Fig. 2, we plot the electronic dispersion 
relation (solid curve) along the symmetry line M-Y-Y of 
the Brillouin zone below the Fermi energy (Ep), when 
si=0.05 and A=l. The Fermi surface of this case is also 



FIG. 3: The positive-negative sign change effect of isotopic 
band shift as a function of the e-ph coupling constant si, on 
a 4x4 square lattice. The ratio AiJk/Ao) at F (upper panel) 
and AI (lower panel) points are shown, where the solid curves 
are calculated by HFA, and the filled circles by QMC. 
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FIG. 4: The changes of lus/l^o (upper panel) and phonon 
spatial variation A(g^) (lower panel) with si, on a 4x4 square 
lattice. The solid curves are calculated by HFA, and the filled 
circles by QMC. 



shown in the inset. By comparing with the results of 
Ref. 16 (dashed curve in Fig. 2), one can see the band 
structure of Bi2212 has been qualitatively depicted by 
our selected parameters. 

In Fig. 3, the HFA (sohd hues) and QMC (filled cir- 
cles) results of the ratio Ai?k/Aw at F [kr = (0,0)] and 
M [kM = (7r,0)] points are presented as functions of si. 
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respectively. Here we see both theories figure out an 
increase of Ai^k/Aw with si, which means if the e-ph 
couphng is strong enough, a small change of phonon en- 
ergy Aw can be amplified to a large band shift A_Ek in 
ARPES. 

In Fig. 4, we show the origin of this anomalously 
large band shift is the phonon softening driven by the 
aforementioned quadratic e-ph coupling. In panel (a), 
the ratio uJs/loq is plotted as a function of the quadratic 
coupling constant Si. Both HFA (solid lines) and QMC 
(filled circles) show that w^/wq declines monotonically 
with the increase of si. Together with this softening of 
phonon energy, the phonon spatial variation A(g^) due to 
the isotope substitution increases gradually [panel (b)]. 
This finally leads to a large band shift, as manifested in 
Eq. (12). 

In Eq. (12) as well as Eq. (27), it is implied that the 
band shift has a momentum dependence, which may re- 
sult in an anisotropy of the isotope effect, and a positive- 
negative sign change of the ratio Aii^k/Aw. As an exam- 
ple, we show this sign reverse in Fig. 3, where panels (a) 
and (b) are for AEi^/Auj at the F and M points, respec- 
tively. Because of the size limitation of QMC, we do not 
have much data for a continuous dispersion along the TM 
direction. Nevertheless, in Fig. 3, both theories clearly 
show that the isotope substitution yields a positive band 
shift at the F point, and a negative one at the M point. 
Thus, one can see the sign reverse of the isotope effect is 
a natural consequence of the momentum dependence of 
the e-ph coupling. 

In Figs. 3 and 4, although at larger si, some deviation 
between the HFA and QMC arises, the overall tendencies 



are in a good agreement, demonstrating that the strong 
off-diagonal quadratic e-ph coupling is responsible for the 
large isotopic band shift. 

Before closing this section, it is worth noting that the 
latest experiment by Douglas et al^ has suggested that 
the unusual isotope effect observed by Gweon et alA 
might be an alignment error of the sample. In contrast 
to the 10-40 meV band shift reported in Ref. 4, new data 
show the shift is only 2±3 meVi^ If this is the case, it 
should be clear that the off-diagonal quadratic e-ph cou- 
pling investigated in this work cannot be very strong for 
the cuprates. 



V. SUMMARY 

In conclusion, by using the HFA and QMC, we study 
the isotope effect of Bi2212 based on an off-diagonal e- 
ph model, where the electrons and phonons are coupled 
in a quadratic form. The HFA enables us to capture the 
essence of present problem, and QMC provides us the ex- 
act results with which the HFA results can be compared. 
Our calculation demonstrates that the isotope-induced 
band shift of ARPES can be well described by the model 
with an off-diagonal quadratic e-ph coupling. We as- 
cribe the large value of A£'k/Atj to the phonon soften- 
ing triggered by the quadratic e-ph coupling. While the 
positive-negative sign change is connected to a momen- 
tum dependence of the coupling. Comparing the results 
of HFA and QMC, we find the agreement is good. 
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